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Preface 


I became interested in the capability of air-to-air missiles first 
hand while on the receiving end of some and the giving end of others 
during my combat tour in Vietnam. When Major Thomas Moriarty 
suggested the subject, I was more than stunntidnatiy interested. The 
need for the development of a maneuvering target came from Dr. Michael 
Caluda, Armament Development Test Center, Eglin AFB, Florida. 

The program was developed using human piiot's observation 
and logic in making decisions. My background is as a tactical fighter 
pilot and instructor pilot. As an instructor pilot, I observed pilots 
make airborne tactical decisions. These observations were incorporated 
in the development of this thesis. 

The assistance provided by my thesis advisor, Major Thomas 
Moriarity, was invaluable. He provided the initial direction and con- 
tinued motivation throughout the development. 

And most of all, I must express my deepest appreciation to my 
wife, Pat, and family for their support and patience during the writing 


of this thesis, without which it is likely I would not have completed the 


project. 
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Abstract 


There was a need for a complete maneuvering target for air- 
to-air combat simulation. A computer program was developed with 
logic and tactical decisions subroutines. The parameters for these 
decisions corresponded to the differences in pilots and relative air- 
craft states. The decision subroutines provided the desired control 
inputs to accomplish the required maneuvers. Validation of the 
simulated target was accomplished by use of five test runs from 


different initial conditions as well as twenty runs from the same 


initial conditions. 


viii 


eo chai mata tt a 


Maneuvering Target Simulation for Testing 
the 


Terminal Guidance of Air-to-Air Missiles 


I. Introduction 


Background 


Successful development of accurate air-to-air missiles depends 
heavily on the results of preproduction simulation testing. For success- 
ful simulation testing and to keep the cost down at all levels of develop- 
ment, a comprehensive and realistic simulation of targets must be 
developed. This will accurately test missile guidance and performance 
before the costly manufacturing of the prototype missile and live-fire 
testing is begun. Many maneuvering targets have been devei ped for 
this purpose and for use in other simulation fields (Ref. 2:51). However, 
a maneuvering target that reacts like a recent graduate of the United 
States Air Force Tactical Fighter Weapons Center (TFWC) to different 
air-to-air missile combat situations is needed to provide an effective 
simulation (Ref. 7). 

The development of the maneuvering target is dependent on the 
aircraft and the pilot. Frequently when maneuvering targets are 
synthesized, the pilot is many times neglected; however, the pilot is 


too important a factor to the results of an air-to-air battle not to be 


included. The outcome of a battle becomes very sensitive to the pilot's 
judgement of range, relative position, and selection of the proper 
maneuver; in fact, these factors sometimes affect the outcome of 
air-to-air combat more than the particular missiles and aircraft 
involved. 

Previous targets have been programmed to follow specific 
tracks or maneuvers for simulated air battles. This may be satis- 
factory for initial testing, but before the costly latter stages of 
testing, a comprehensive air-to-air combat simulation must be 


available. | 


Objective 
The development of a fully maneuvering target was the 
objective of this thesis. The requirements of the maneuvering target 
are: 
1. Respond with the desired maneuvers to the attack of the 
missile. 
: 2. Include a human pilot range-judgement error. 
3. Update the desired maneuver and the inputs throughout 


the simulation, 


Description and Scope 


The maneuvers incorporated in the simulation should be the 
latest in evasive tactics. The program must be able to accommodate 
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additional or new maneuvers with a minimum of programming changes. 

The development of the maneuvering target simulation resulted 
in the program not being limited to any specific aircraft or missile. 
The program can be easily adapted to accommodate most aircraft and 
missiles. 

The maneuvering target uses decision logic to respond to the 
given aircraft-missile situation, to determine the finite control inputs 
to fly the evasive maneuvers, and to continuously update the decision 
process so as to select an evasive response. 

A realistic maneuvering target must not be optimized because 
a pilot has only one chance at each combat situation. During air-to- 
air combat and especially in maneuvers against launched missiles, 
the pilot has a very limited time to estimate the range and closure 


rate and then choose a maneuver. Asa result, there are incidents 


a 


where the pilot selects other than the optimal maneuver or even the 


wrong Maneuver, 


The responsive target was developed in the following manner: | 


Dink as a 


1. The maneuvers were selected for the program. 
2. The basic decision logic was developed. 


3. The maneuver noise was selected. This was the marner 
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in which the pilot's response was incorporated into the simulation. 


4. Control inputs were developed for the maneuvers in 


step 1. 


tt th 


5. The simulation was then programmed and validated. It 
was necessary to run tests from different positions as well as 


numerous runs from the same initial conditions. 


Organization 

The thesis includes the combat scenario and description, 
analysis of the computer program, validation of the aircraft model, 
plus results and conclusions. The program listing, printed results, 
data input, development of the ae SRE coefficients, and the 
missile model are included in the appendices. 

The combat scenario, Section II, includes the development 
of the mathematics of the aircraft model. In addition, the develop- 
ment of the decision logic and missile seeker noise are presented. 

Each individual function and subroutine are discussed in 
Section III, Computer Program. The calling sequence of the rouiines 


and the theory behind the development of the main program and sul- 


programs is covered, 


RST ener... 


In Section IV, the aircraft model is validated as well as the 


decision logic. The initial conditions, aircraft characteristics 


Pe hep PETER 


and the equations of motion are developed. 


The results of the five validation runs and the 20 runs from 


the same initial conditions are in Section V, Results and Conclusions. 


Each validation run is considered individually as well as the total 


results of the twenty consecutive runs. 


>» 


Il. Combat Scenario and Description 
pa ba lathe ean cbse Neteen cesta dated wdlglah ahs 


Computer simulation of air-to-air combat is a means of pro- 
viding for testing and evaluating tactics, maneuvers, and missile 
hardware. When the simulation is designed to follow pre-selected 
maneuvers or tracks, the operator in reality dictates the results. 
The best maneuver or tactic is a result of trial and error by the 
simulation operator. 

The optimal type of simulation is not realistic. Each situation 
must allow only one selection by the pilot. In addition, maneuvering 
noise must be introduced for errors in judgement of the pilot. This 
idea will be developed further in the discussion of Subroutine Pilot. 

Combat logic was developed from the relative states of the 
two vehicles and this logic provides feedback to the decision making 
process. The feedback continuously updates the decision parameters 
so that the maneuvering target will react as if controlled by a human 


pilot. 


Model Description 


The aircraft is represented by a center of mass model. The 


aircraft contains a body centered coordinate system, with the origin 


located at the center of gravity, C.G., of the aircraft. The Xp axis 


is directed out the front of the aircra/ Y¥37axis is directed out 


the right wing of the aircraft, and the Zp-axis completes the right handed 
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orthogonal coordinate system, therefore, is directed out the bottom 
of the aircraft. The orientation of the body axis is shown in Figure 1. 

The missile is described by the same type of body axis 
coordinate system. The pilot superimposes a spherical coordinate 
system over the body axis of the aircraft so as to relate the missile 
to his own frame. This relationship is shown in Figure 2. When 
range (R) is utilized, it refers to the distance between the aircraft 
and the missile, Zeta (¢) is the asimuth angle to the missile from © 
the aircraft. The angle is measured positive right in the Xp-Y, 
plane of the aircraft. Eta (1) is the elevation angle from the aircraft 
to the missile with positive being up from and perpendicular to the 
Xp-Yp plane (Ref. 1:9). 

The navigation frame has its origin at the C.G. of the aircraft. 
The Yy-axis is directed toward east and the Zyn-axis is directed 
downward along the gravity vector, g. The Xy-axis completes the 
right hand orthogonal coordinate system. The Xyraxis is directed 
north. 

The wind axis reference frame is also used. The origin of 
this frame is the C.G. of the aircraft. The Xyraxis is directed 
along the velocity vector of the aircraft with respect to the atmosphere. 
The Zy7 axis is in the plane of symmetry of the aircraft and is 
directed perpendicular to the Xyr axis. The Yyw axis completes the 


right-hand orthogonal coordinate system (Ref. 2:109). 
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Before describing the final coordinate system that is utilized 
in the aircraft simulation, a few basic assumptions have to be made: 

1. The earth's fixed frame is considered an inertial frame. 

2. The earth's rotation is neglected. - 

3. The earth is assumed to be flat. 

4. ‘Greedy is assumed to be constant in direction and 
magnitude. 

5. Air mass movement is assumed to be zero. 

With these assumptions included, the earth's surface fixed 
coordinate frame is introduced. This frame parallels the aircraft 
navigation fieines The origin of the earth's surface frame is fixed 


on the earth's surface and is directly below the aircraft C.G. at 


time zero. At the beginning of each engagement, the aircraft's 
horizontal coordinates are zero and the vertical coordinate is equal 
to the negative altitude of the aircraft. It can be seen that at time 
zero any vector in the navigational frame will have the same orienta- 


tion in the earth's surface reference frame and vice versa (Ref. 1:10). 


Decision Logic 


Maneuver selection is based on the instantaneous relative 
states between the aircraft and missile. The decision logic con- 
sidered the range between the two vehicles and the closure geometry. 


To transform the decision logic into the reactions of a human, the 
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selection parameters used by the pilot are incorporated into the decision 
logic routine. A pilot will see the missile as two angles and estimated 
range. The pilot describes the missile as, "Missile, five o'clock high 
at 5000 feet. '' Converted into mathematical terms, the position trans- 
lates as range of 5000 feet, azimuth 150 degrees, and elevation of 

45 degrees (Ref. 1:10). 

At long ranges, the assumption of point mass is selected due to 
the inability of the pilot to determine relative sizes. He can estimate 
range, azimuth, and elevation. 

After determining an estimated position, the maneuver selection 
is based on the relative state vectors. The state vectors are resolved 
into two two-dimensional angles. First, the aircraft angle-off is 
considered. The angle-off (9 off) is defined as the angle between the 
line-of-sight vector and the missile's velocity vector. The cone angle 
(8 cone) is the angle between the velocity vector and the line sight - of - 
sight vector of the aircraft. 

Over a period of time, the pilot is able to estimate the rate of 
change in the relative position between the aircraft and the missile. 
Therefore, the pilot adds the estimated range rate (R) azimuth rate (¢) 
and elevation rate, tl. to his knowledge for use in maneuver selection 
(Ref. 1:11-12). 

A pilot is now able to make a logical maneuver selection based 


on; 
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R - Range R - Range Rate 

€ - Azimuth ¢ - Azimuth Rate 
7], - Elevation 1 - Elevation Rate 
~ - Angle Off @ - Cone Angle 


The relationship between these variables are shown in 
Figure 3. 

Before the previous eight parameters can be determined, the 
earth's fixed coordinates have 0 be transformed into the aircraft 
body fixed frame. The aircraft states consisted of the following 
components: 

| X, Y, Z - position 
u, v, w- velocity 
Y, 8, » - Euler angles (heading, flight path angle, 
and bank angle) 
The missle states consisted of the same components. 
The missile's relative position is converted in the aircraft's 


navigation frame. The velocity components are also converted. 


XMF" XM 7 *P ‘1) 
| YMF = Yq ~ Yr (2) 
7MF~ 7M 7F (3) 
Uye = Uy Up (4) 
Yur = Ym ~ ‘pr (5) 
w = wWlo-w (6) 
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The missile position and velocity components in the aircraft's 
navigation frame are converted into the aircraft's body fixed frame by 
use of the Euler angle transformation. 

c8c¥ césY -s@ 


TB = |s¥s@c¥-cysY sosés¥t+cys¥ svcd (7) 
cpsOcl¥+sys¥ cPs&s¥-sgc¥ cypce 
~— 


where gs denotes sine and c denotes cosine. 


Xu | *ME a 

Yup . TB Yur (8) 
2MB “MF 

“MB MEP | 

YMB “Typ | YMF (9) 
YMB | | 


The azimuth and elevation of the missile with respect to the 


aircraft are determined in the following manner: 


si 
Smr~ 5" | Yup (10) 
R 
Tare = sine? [-z (11) 
MEF F 
4 


The simplest means of determining the range between the 
missile and aircraft is to use the position coordinates of each in the 
earth's surface fixed frame, The range is computed as the square 
root of the sum of the squares of the differences of the three 


components of position. 
14 
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R = X + ¥ sled ot (12) 
MF MF ME 


The closure rate or range rate is computed by use of differ- 
ences in the relative positions and the differences in the components 
of velocity. Both velocity and position components are in the air- 
craft body fixedframe. By use of ei eetetnst substitution, the 
range rate is determined by use of the velocity differences and 
two angles, azimuth and elevation. The range rate is computed as 


the difference of the missile and aircraft velocity components along 


the line-of-sight. 


: + Z ; 
R* "v5 *mp*’up’mp*”’mp’mp (3) 
R 
or 
R= Uy COS Tare cos CMR + VB cos Ter sin Cer 
Ss: Wp sin Tae (14) 


The azimuth rate and elevation rate are then determined. To 
determine these variables, the heading angle of the missile, expressed 


in the aircraft's body fixed frame, ve) is derived first. Yar is the 


angle between the aircraft's x axis and the velocity vector component 


of the missile in the XBY,, plane of the aircraft, 


u 
‘my * cos! —y-—ME—, 


(vet Yup) 
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If vyp iS negative, MEF becomes the negative of the above 


equation. 


(16) 
2 2 ! 
[Ue * Yup i? 


Another angle of interest is the flight path angle of the missil« 


in the X,Y jp plane of the aircraft, @\i 


smd -WMB 
6... = Sin a (17) 
MF 2 + y2 Z 1 


+ : 
(“mp ’msp. ” MB]? 


With YMEF and ouE defined, the azimuth rate €, and elevation 


rate, Nygpr are then determined. 


2 De ok 
1 “an ny 
Se taal 2280/2 * Gee “Tage 
MF R cos Nya 
: 2 2 
IMF [ump + VB jz cos (Ye - Cure sin Nag Yypp ©°S 1 


R 
Refer to Figures 4 and 5 for the orientation of C and Th 


The next angles to be determined are the angle off, Pore, and 


cone angle, %.one- The scaler dot product of R and the aircraft's 


velocity vector, V, determine the Pore: The mathematical solution: 


cos (orp) = V_. R (20) 
Vv R 


vi 
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Mir (19) 


cos (Poff) = MF 444 + YM YuMr’m + ZK Wu (21) 
VR 

The angle-off is defined as the angle between the line-of- 
sight vector and the missile's velocity vector. 

The cone angle is defined as the angle between the aircraft's 
velocity vector and its line-of-sight vector. The mathematical 
value is determined by: 

cos (8 ne? = cos Cmr cos Vr (22) 


(Ref. 4:77-83) 


Missile Seeker Noise 

The calculation of the desired variables for the missile are 
the same as for the aircraft except the relative positions are 
reversed. Another area must be introduced at this point. The 
missile seeker noise is not just a constant Gaussian noise source. 
A scaling factor is introduced to increase or decrease the standard 
deviation, 6, of the guidance error of the missile as a factor of the 
missile's and aircraft's relative orientation. 

First, the radar guided missile noise is assumed to be 
minimized when the missile is directly in front of or behind its 
target. As the relative line-of-sight rotates around to the side of 
the aircraft, the azimuth guidance error is assumed to be the 


largest. The same phenomenon occurs in the elevation errors. The 
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point of the largest guidance error due to elevation is assumed to be 
when the missile is directly above or below the aircraft. 

The seeker noise for the infra-red guided missile is modeled 
in the same way except for the position of maximum guidance error. 
The position of maximum noise is assumed to be directly in front of 
the missile's target, while the minimum noise is assumed to be when 
the missile is directly behind the target. The infra-red guided 
missile tracks the maximum heat point of the plane following the 
aircraft. As the relative position rotates from directly in front of 
the target to the rear, the scale factor decreases, therefore, the 
guidance error decreases. 

It was assumed for this simulation that the missile tracks 
the plume of the aircraft. Therefore, an additional X-axis guidance 
error is present and it contributes significantly to the X-guidance 
error. 

The angles used to determine relative position are developed 
in the description for Subroutine ATTACKF. These angles represent 
the relationship of the missile to the aircraft. The angles deter- 
mine the actual presentation the missile seeker attempts to track. 

The X-axis guidance error is assumed not to be entirely a 
result of the relative azimuth angle. As the relative elevation angle 
increases, the X-axis guidance error also increases. The same 


is true of the Y and Z guidance errors. Therefore, it is assumed 


20 


ee et NELLA LALOR ALL ETO ta ane ee OO em 


that elevation and azimuth relative angles effect all three of the guidance 


errors. 


To define the scale factor as a non-dimensional variable, the 


relative angles are divided by the inaximum value that the angle could 


attain. The scale factor is used to compute the standard deviation of 


the Gaussian distribution. The Gaussian number generator routine is 


taken from Gordon (Ref. 3:114-115). The actual guidance error is 


computed so that the deviation from the previously computed error is 
small. The noise is a high frequency noise that is integrated out by 
the missile and, therefore, has little effect on the missile guidance. 


The actual modeling of the missile and associated noise is beyond the 


scope of this thesis. 
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III. Computer Program 


Program Evasion 


The computer program was not only developed for the accuracy 
of the simulation but also simplicity, size, and modular construction. 
This led to reasonable memory requirements, compilation and execu- 
tion times. 

Arrays are used for initial conditions, aircraft and missile 
data, and all variables derived in each subroutine. This provides for 
ease of transfer between the main program and subroutines and keeps 
the program variables simple. 

The versatility of the program is extended by use of an integer 
to denote non-responsive or responsive targets and infra-red (IR) or 
radar guided missiles. Other type guidance systems could be added 
without a great deal of change to the basic program. 

The main program is a series of calls for subroutines. For 
each run, the initial conditions and aircraft/missile data routines are 
called once. Subroutines are called once at the beginning to change 
all angles/angular rates from degrees to radians. 

The subroutines called during each integration step will be 
discussed separately. In general, these subroutines determine and 
update the parameters, select the desired maneuvers and control 
inputs and print the desired outputs at regular intervals. At the com- 


pletion of each integration step, conditions are compared with 
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maximum values. A flow chart is depicted in Figure 6. 


Data Inputs 

Subroutines DATAFTR, DATAMIS, INITFT, and INITMIS read 
in the respective aircraft and missile data, along with the initial 
conditions of both vehicles. These four subroutines are called once at 
the beginning of each run by Program EVASION. 
Subroutine ANGLESF/Subroutine ANGLESM 

The angle and angular rate conversion subroutines for both 
vehicles transform the degree inputs into radians by use of the 
Function DEGZRAD. ANGLESF and ANGLESM are called following 
the data inputs in EVASION. 

Subroutine PRINT 

Subroutine PRINT performs three tasks. It calls Subroutines 
ATTACKF and ATTACKM for the respective vehicles, converts 
desired output variables from radians to degrees, and prints the desired 
output variables in the proper format (Ref. 1:73). 

The printed output provides for easy comparison between the 
position, velocities, and orientation of the two vehicles plus the 
relative position of the missile with respect to the aircraft. 

Subroutine PRINT converts the desired program outputs from 
radians to degrees by use of the Function RADZDEG. 

Subroutine PRINT calls the Subroutines ATTACKF and ATTACKM 


for the respective vehicles. These subroutines will be explained later. 
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' Subroutine PRINT is called at three points in the simulation. 

It is called at time zero to initialize the relative position variables 
and print out initial conditions. PRINT is called at the completion of 
each integration, so that the desired outputs will be printed at regular 
intervals throughout the simulation. At termination, PRINT is called 
to record the end conditions. 
Subroutine ATTACKF/Subroutine ATTACKM 

Subroutines ATTACKF and ATTACKM determine the relative 
position and velocity components of the opposing vehicle in its own 
navigational frame. These components are then transformed into the 
respective vehicle's body fixed frame. 

The relative angles are computed during each integration upon 
which the desired maneuvers are selected. 

Subroutines ATTACKF and ATTACKM are called in Subroutine 
PRINT. ATTACKM calls the two seeker noise silt: IRNOS 


and RADNOS, which will be discussed later. 


Subroutine PILOT 


The maneuvering target must react like a human pilot. It will, 
therefore, be less than perfect in its judgement and estimation of 
range. By use of a random number generator, a range judgement 
error is introduced. The judgement error allows for variation in the 
pilot's ability to estimate range. 

There are two decision points in the simulation; 6000 and 1000 


feet are the ranges at which new maneuvers are selected. The range 
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errors are bounded by specific values which are determined from 


flying experience. If better values are available, they should be in- 
cluded. At the 6000 feet point, the range errors have a uniform distri- 
bution = 1000 feet, while at the 1000 feet point, the distribution is 
within + 500 feet. In turn, this range judgement error changes the 
execution point of the desired maneuvers. The ranges are then used 
as parameters in Subroutine EVAMANU to select the desired maneuver. 

This subroutine is an over-simplification of the judgement of a 
human pilot but it does provide for a variation in pilot skill. All levels 
of skill are present in every flying unit. 

PILOT is called at the beginning of each run in Program 
EVASION. 
Subroutine EVAMANU 

Maneuvers against a launched air-to-air missile must be 
simple in execution. Generally, the defense against a missile fired 
in the aft quarter is a maximum turn rate "break" turn into the plane 
of the missile. As the missile closes to approximately 1000 feet, a 
rapid roll out of the plane is initiated. Obviously, because of closure 
rates, the timing of this last portion is very subjective. 

The attempt is to initially generate maximum line-of-sight rate 
for the missile and then to further complicate plane corrections during 
end game maneuvering. This defense works well against rear hemis- 


phere missiles launched near the heart of the firing envelope. 
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Long range firings provide another problem. A hard turn into 
the missile only supplies closure so that the missile is not required to 
fly a pursuit type curve. Therefore, a hard turning target is hit from 
the front quarter rather than the tail. 

There are two options available against long-range missile 
firings: 

1. Hard turn away to put the aircraft outside maximum range. 

2. Turn away enough to force the missile into a "tail chase" 
and then break hard into it at 6000 feet, hopefully causing an overshoot. 

Head-on or beam shots are best handled by a hard pulling turn 
to put the aircraft outside of range. However, if the pilot does not 
want to lose sight of the missile, a hard turn away should be initiated 
to force the missile into a pursuit type curve and proceed with an 
overshoot. 

At any time, breaking downward at high calibrated airspeeds 
to lower altitudes not only shrinks the missile envelope, but also 
increases background noise and clutter. 

Early turns into radar guided missiles do not help. These 
turns just present a larger radar return and do not affect the 
maneuvering required by the missile. 

Near the maximum firing range, the aircraft should extend 
away from infra-red guided missiles. This is extremely hard to 


judge and the normal response for IR missiles is to turn into the 
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missile and increase angle off, orf, the tail to either put it outside 
source limits or to increase tracking and maneuvering difficulty. 

Infra-red guided missiles are very fast and maneuverable. 
The best way to counter launches from within the heart of the firing 
envelope is to have the break be in the plane of the attacking aircraft's 
wings at launch. This turn can be increased to a break after launch, 
making the initial turn in the plane of the missile. 

Subroutine EVAMANU is called in Program EVASION normally 
one to three times depending on the initial range between missile and 
aircraft. EVAMANU is called at the 6000 and 1000 feet decision 
points as well as at the initial point of the simulation. 

Subroutine THRUSTF /Subroutine THRUSTM 

Subroutines THRUSTF and THRUSTM normally compute the 
thrust of the respective vehicles by use of Subroutines ATMOS and 
TBLNOC (Ref. 5). Maximum, Military, and Idle thrust are deter- 
mined as well as the present thrust of the aircraft. In both thrust 
subroutines, the total velocity, mach number, and dynamic pressure 
of each vehicle are determined. The missile is considered a constant 
thrust and velocity vehicle for the test runs of the simulation. For 
actual missile testing, THRUSTM would be used to compute the 
thrust of the missile. 

THRUSTF and THRUSTM are called in Program EVASION 


during each integration. 
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Subroutine DESINP 

Subroutine DESINP provides the desired control] inputs for the 
selected maneuvers from EVAMANU. ‘The inputs are updated during 
each integration. The control inputs are angle of attack, sideslip 


angle, angie of bank, and desired thrusi. These control inputs corre- 


spond to elevator, rudder, aileron, and throttle control (Ref. 4:103- 


104). 

Program EVASION calls Subroutine DESINP during each 
integration. 
Subroutine PURSUIT * 


Subroutine PURSUIT provides the control inputs for the missile 
in the same manner as DESINP did for the aircraft. PURSUIT provides 
the desired heading and pitch for the missile. 

If, at any time, the missile goes ballistic, PURSUIT includes 
a test condition to return to the main program. Therefore, the desired 
inputs remain the same. 

PURSUIT is called during each integration by Program EVASION. 
Subroutine INPUTSF/Subroutine INPU'SM 

Finite control rates are supplied to the simulation by Subroutines 
INPUTSF and INPUTSM to their respective vehicles. These instantaneous 
contro] rates are required to make the simulation realistic. The control 
rates are functions of the control errors. In INPUTSI’, angle of attack, 


back angle, sideslip, and thrust are computed. The pitch and heading 
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angles of the missile are determined in INPUTSM. 

/ Using angle of attack (@) as an example, the error is deter- 
mined between the desired angle of attack (a) and the actual a of the 
aircraft. 


ERROR @ = apy @ (23) 


The actual control inputs are determined for the simulation by use of 


the following limits: 


ERRORa 4@« = Py OMAX feo w \ 


it SP Max (24) 
(At) Po Ora x(t) 


where: 

@ - rate of change of angle of attack 

Py, 9 ~ guidance parameters se aircratt 

Or AX - max rate of change of angle of attack 

At - integration interval 

Atte’ & was determined, the new qa for the aircraft is com- 
puted: 

® = Ay! 1 = (SIGN ERROR Q) @ (At) (25) 
The relation of control rates to control error are illustrated in 


Figure 7 (Ref. 4:93). 


Again using angle of attack as an example, the relationship becomes: 


a= wy) (Ref. 4:94). 
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The finite control inputs can be by-passed if the error is zero. 
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INPUTSF and INPUTSM are daltad in Program EVASION during 
each integration. INPUTSM is also called in FORCESM if either the 
vertical or hotizontal g forces exceed the maximum for the missile. 
This call provides for a recomputing of the pitch and heading angles. 
Reduced angular rates are supplied by FORCESM which maintain the 


g forces within limits. 


Subroutine FORCESF/ Subroutine FORCESM 


Subroutines FORCESF and FORCESM determine the forces and 
moments acting on the respective vehicles. By use of Function TBLNDC, 
the coefficients and derivatives are determined in FORCESF for the air- 
craft. The total coefficients were computed and the total forces and 
moments are then determined. 

Other parameters defined in FORCESF are the pitch, bank, and 
sideslip angles in the wind coordinate cians 

In FORCESM, only the vertical and horizontal g forces are com- 
puted. This is due to the very simple missile supplied to the simulation. 
In actual tests, FORCESM should be expanded to the same extent as 
FORCESF. 

FORCESF and FORCESM are called during each integration in 
Program EVASION. FORCESM calls INPUTSM if either the horizontal 


or vertical g forces are exceeded. There is fur:.er discussion in 


Appendix D and Section INPUTSM 
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Subroutine UPDATEF/Subroutine UPDATEM 

The actual integration is accomplished by use of Subroutines 
RKDESF and RKDESM. The deferential equations are supplied to the 
RKDESF and RKDESM routines by subroutines F and M respectively. 

UPDATEF and UPDATEM calls RKDESF and RKDESM respec- 
tively during each integration. Program EVASION calls UPDATEF and 
UPDATEM during each integration. 
Subroutine IRNOS/Subroutine RADNOS 

Subroutines IRNOS and RADNOS are routines that determined 
the guidance error for infra-red and radar guided missiles. The basic 
Gaussian random generator is modified by a scaling factor that changes 
the dispersion of the noise as a function of the relative orientation 
between the missile and the aircraft. The scaling factor is combined 
with the standard deviation to produce a corrected value as a function 
of the relative orientation. The missile aim point deviation is computed 
in the standard manner. Further discussion of the methodology in the 
noise routines is in the section ''Missile Seeker Noise" in Chapter II. 

IRNOS and RADNOS are called by ATTACKM according to the 


type missile being simulated. It is called during each integration. 


Subroutine REDUCE 


Subroutine REDUCE is a simple program to reduce the time 
rate of change of the range between the two vehicles. When the range 
covered during one integration interval is equal to or less than the range 


between the two vehicles, the integration time interval is reduced so 
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that the range covered in one integration interval is less than one foot. 
This procedure is a simple means of providing accuracy of less than 
one foot during the termination phase. REDUCE is called in the Program 
EVASION. 
Miscellaneous Subroutines and Function Subprograms 

The previously discussed programs and subroutines are the main 
solution to the simulation problem. There were several functions and 
subroutines developed to perform specific calculations or transforma- 
tions. In addition, two subroutines and a function subprogram are used 
from the Air Force Institute of Technology (AFIT) Subroutines Library 
(Ref. 5). 
Function TBLNOC 

This function is from the AFIT Subroutines Library. It is an 
n-dimensional table look-up function that is used in THRUSTF and 
FORCESF for aerodynamic coefficients and thrust parameters. Aero- 
dynamic coefficients are stored as a function of angle of attack and 
Mach Number. Thrust parameters are stored as a function of Mach 
number. TBLNDC is from the AFIT Subroutines Library (Ref. 5) and 
is furnished in Appendix A. 
Subroutines RKDESF and RKDESM 

These subroutines are taken from the AFIT Subroutines Library. 
RKDESF and RKDESM are fourth order Runga-Kutta differential 


equation solvers for the respective vehicle. The subroutines RKDESF 
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and RKDESM are called in the respective UPDATEF and UPDATEM 
subroutines (Ref. 5). and is furnished in Appendix A. The differential 
equations are supplied by subroutines F and M for the respective inte- 
gration routines. 
Subroutine ATMOS 

This subroutine is taken from the AFIT Subroutines Library. 
ATMOS is a tabulation of the 1959 ARDC atmospheric tables. The 
atmospheric parameters are determined as a function of altitude by 
linear interpolation. From ATMOS, the speed of sound, the density, 
and the density ratio are used to calculate many variables in the 
simulation. ATMOS is called in THRUSTF, THRUSTM and INITMIS 
{Ref. 5) and is supplied in Appendix A. 
Subroutines F and M 

These subroutines are called by RKDESF and RKDESM respec- 
tively. They contain the differential equations of motion to be solved 
by RKDESF and RKDESM. 
Subroutine TRANN2B 

Euler angle transformation equations are used in TRANN2B to 
change a vector from the vehicle's navigation frame into the vehicle's 
body frame (Ref. 2:116-117). 
Subroutine TRANB2W 

By use of the angle of attack and sideslip angles, TRANB2W 


transforms a vector from the vehicle's body fixed frame into the 
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vehicle's wind fixed frame (Ref. 2:116-117). 


Function DEG2RAD 

The angles and angular rates are normally inputed as degrees. 
DEG2RAD transforms degrees into radians. 
Function RAD2DEG 

For convenience, the angles/angular rates are printed in 
degrees. After computations, RAD2DEG transforms radians into 
degrees. 
Others 

The general utility functions are taken from the FORTRAN 
Extended Library (Ref. 6:8.1-8.12). They include: 

1) Intrinsic Functions - IFIX, SIGN AMAX 1, AMIN 1, ABS, 

MOD 
2) Basic External Function - SIN, COS, SQRT, ACOS, ASIN 


3) Utility Subprograms - EOF, RANP 
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IV. Aircraft Model 


Validation 

Validation of the aircraft simulation model was accomplished 
in basically two ways. The flying experience of the author was used 
as well as the Technical Order of the F-4C (Ref. 8). 

After the complete aircraft model was developed, the simula- 
tion was accomplished with a missile with zero closure rate to 
validate all of the programmed aircraft maneuvers. Using the author's 
flying experience, the performance of the aircraft was determined to 
be consistent with present day fighter aircraft. 

The simulated performance was compared to published per- 
formance data in the F-4C Technical Order (Ref. 8:A9-94, 95). 
Performance parameters used were rate of turn, angle of bank, and 
radius of turn as a function of altitude and Mach number. Comparison 
was very favorable, Deviations were a result of simplification of the 
total aeiedveanie coefficient computations in FORCESF. For com- 
parison, Table I displays some test and technical order results. 

The decision logic of the aircraft required validation. During 
the zero closure rate missile simulations, the maneuver selection and 
aircraft reaction and maneuvering were compared to tactical flying 
manuals and the author's combat flying experience. The aircraft per- 


formed as if controlled by an average fighter pilot. 
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Initial Conditions 


The initial conditions of the aircraft consist of: 
Position coordinates (xp, Ye zy) 
2. Velocity components (Uy» Vie Wy) 
Orientation angles (heading, bank, flight path) 
Angular velocity components (p, q, r) 

p - about XB axis 

q - about Yp axis 

r - about 23 axis 
The initial conditions 1 and 2 are in the earth surface fixed 


reference frame. Number 3 is in the aircraft body fixed reference 


frame. 


The initial conditions are used in the first integration for 


maneuver selection and computations. 


Aircraft Characteristics 


The characteristics of the aircraft are required to determine 


the coefficients for the equations of motions. The aircraft chosen to 


“se nea RE 


be simulated was basically the F-4C. The actual aerodynamic coefficient 


aon eee 


equations and stability derivatives used in FORCESM are in Appendix C. 


The aerodynamic coefficients consist of the basic coefficients and 


stability derivatives. They were compiled as functions of the Mach 


Number and angle of attack. The thrust parameters were compiled as 


a function of Mach Number alone (Ref. 9). 
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Equations of Motion 


The differential equations of motions used in the simulation are 
taken from ETKIN (Ref. 2:149-150). These equations were derived with 
the following assumptions: 

1. Rigid body 

2. Piane of symmetry in combined wind and body axes. 

3. Flat earth 

4. Constant mass 

5. Constant gravity acceleration 

“ 

6. Atmosphere is at rest relative to the earth. 

7. Centripetal acceleratiofi associated with earth rotation is 

neglected. 


The differential equations of motion of the maneuvering target 


J 
are: 
al Se 
v= m(Tyw - D - mg sin 6 w) (26) 
5: Regd 
P = Iyy [L - Iyz (i + pq) + (lyy- Ip 7) qr) (27) 
pie 2 2 
ces 
8w = qw cos gw - rw sin gw (30) 
Yw = (qw sin gw + rw cos ww) sec Ow (31) 
| Xp = V cos @w cos ¥w (32) 
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) V cos ®@w sin 4 w (33) 


Zp = -V sin ¥ w (34) 


The equations for @w, a, 8, (Ref. 2:150) are not used in Sub- 
routine F. The finite control inputs ow, qa, 6, are developed in Sub- 
routine INPUTSF. 

Integration of equation 26 produces the force equation in the x 
direction. The force equations for the y and z directions are not 


differential equations. Forces in the y and z direction are computed 


in FORCESF. 
Fy = Tw - Sideforce - mg cos Sa sin Py (35) 
zo * Tow - Lift + mg cos @y, cos Pw (36) 


V Results 


Five test combat battles were simulated from different initial 
conditions of the missile, different airspeeds of the aircraft, different 
types of missiles and aircraft. In addition, production runs were 
accomplished by repeating the same battle twenty times from the same 
initial conditions, type missile and aircraft. 

The aircraft started at the same point in each battle; x = 0, 

y = 0, z= -20000. The velocity of the aircraft increased by 100 feet 
per second for each battle. The first aircraft had an initial velocity of 
700 feet per second. All angles and angular rates are zero including 
heading which was north. 


The two types of missiles, infra-red and radar guided, were 


ote 


simulated. Each type was directed once against a non-responsive air- 
craft for initial testing. The other three battles were from three dif- 
f ferent ppsitions to further test the aircraft simulation. 

The initial conditions of the missile and aircraft are listed in 
Table Il. Aircraft 3 and 4 were used with missile 3, which was posi- 
tioned in front and behind respectively, for the production runs. Actual 
test numbers and results are listed in Appendix B. 
Battle 1 and 2 

The first two missiles were flown against non-responsive air- 
craft. Both type missiles scored hits on the target. The guidance 


errors were small due to the small angle off in azimuth and zero 
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elevation angle. During actual missile testing this type of run would 
be used to calibrate and initially test missile systems and simulation. 
Battle 3 

The radar guided missile was positioned for a 30 degrees off 
forward pass. The aircraft was a responsive target. Because the 
range was greater than 6000 feet and the missile was in front of the 
aircraft,selected a SPLIT S as the first maneuver. 

At the 6000 feet range point, EVAMANU updated the maneuver 
selection to a HARD BREAK so as to force the missile into a aft 
quarter attack. The turn was then reversed to atterapt to force the 
missile into an overshoot. 

A VERTICLE DIVE FOLLOWED BY A HARD PULL UP was 
selected at the 1000 feet point by EVAMANU. The aircraft was 
successful in evading the missile. The guidance was again minimal 
due to zero elevation angle and a small angle off in azimuth. The 


results are shown in Figure 8. 


Battle 4 


The infra-red guided missile was flown against a responsive 
target. The missile was positioned below, behind, and to the left of 
the aircraft. For this situation, EVAMANU selected a VERTICLE 


DIVE away from the missile due to the range being greater than 6000 


feet. 
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At the 6000 feet decision point, EVAMANU changed the maneuver 
to a HARD BREAK into the missile. The aircraft attempted to cause an 
overshoot by trying a SPLIT S at the 1000 feet range point. 

The aircraft was successful in evading the missile. The guidance 
error was relatively large due to large elevation and azimuth angles. 
The results are plotted in Figure 9. 
Battle 5 

The radar guided missile was directed against a responsive 
target. The range was greater than 6000 feet. The missile was below, 
behind, and to the right of the aircraft. EVAMANU selected a VERTICLE 
DIVE away from the missile. 

EVAMANU updated the selected maneuver at the 6000 feet point 
to a HARD BREAK into the missile. The break continued until over- 
shoot or 1000 feet range, whichever came first. EVAMANU selected 
a SPLIT S at the 1000 feet range point. The maneuver was continued 
until the missile scored a hit. As in battle 4, the guidance error was 
relatively large for the same reasons. The results are shown in 


Figure 10. 


Statistical Characteristics 
Twenty consecutive runs were made from the same initial condi- 
tions to test the stocastic properties of the guidance errors and pilot 


judgement errors developed in the program. The initial conditions are 
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from the aircraft in tests three and four. The missile's initial conditions 
are those of missiles from test battles one and three. The conditions 
are listed in Table II. 

Using the infra-red guided missile from the battle one and air- 
craft from test four, the resulting statistics are recorded in Table III. 


The results of the missile and aircraft of battle three are recorded in 


Table IV. 
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Table III-1 


2 


The statistics of the following parameter are the results of twenty 
production runs from the same initial conditions. 


Decision Point 1 Range (feet): 
Mean= 5999. 90 Standard Deviation= 629. 20 
Minimum= 5025. 00 Maximum= 6950. 00 


Median= 6210.00 


Decision Point 1 Time (Seconds): 


Mean= . 90 Standard Deviation= e2e 
Minimum= ~47 Maximum= 1.35 
Median= . 81 
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Table III-2 
The statistics of the following parameter are the results of twenty 
production runs from the same initial conditions. 
Decision Point 2 Range (feet): 
Mean= 1150. 90 Standard Deviation= 274. 39 
Minimum= 592.00 Ma ximum= 1453. 00 


Median= 1284.50 


Decision Point 2 Time (seconds): 


Mean= 3.14 Standard Deviation= -15 
Minimum= 2.98 Maximum= 3.50 
Median= 3.10 
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Table III-3 


The statistics of the following parameter are the results of twenty 


production runs from the same initial conditions. 
Time to go from decision point 1(seconds) 

Mean= 2.74 Standard Deviation= 
Minimum= 2. 31 Maximum= 


Median= 2. 85 


Time to go from decision point 2(seconds) 


Mean= - 50 Standard Deviation= 


Minimum= .16 Maximum= 


Median= «59 
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Table Il-4 


Termination Time Time To Go Time To Go | 

From Point 1 From Point 2 | 

(seconds) (seconds) (seconds) | 
3.61 3.04 . 63 
3.64 2. 94 wet 
3.65 2. 85 .38 
3.60 3.13 .59 
3.66 2.33 16 
3.65 2. 84 . 63 
3.63 2. 98 . 64 
3.66 2. 36 -53 
3.66 2.31 . 66 
3.64 2.92 .65 
3.61 3.09 37 
3.67 2.65 .31 
3.66 2.34 67 
3.63 2.94 -61 
3.67 2.56 47 
3.67 2.54 . 60 
3.64 2. 89 57 
3.66 2.50 42 
3.67 2.65 ~45 
3.62 3.00 - 40 
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Table IV-1 
The statistics of the following parameter are the results of twenty 
production runs from the same initial conditions. 
Decision Point 1 Range (feet): 
Mean= 5952. 70 Standard Deviation= 668.68 
Minimum= 5011.00 Maximum= 6956. 00 


Median= 5700. 00 


Decision Point 1 Time (seconds): 


Mean= £96 Standard Deviation= Pom lay 
, 
Minimum= 1. 30 Maximum= 1. 81 
Median= 1.63 
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Table IV-2 


The statistics of the following parameter are the results of twenty 
production runs from the same initial conditions. 


Decision Point 2 Range (feet): 


Mean= 937. 30 Standard Deviation= 297.62 


Minimum= 500. 00 Maximum= 1427.00 


Median= 963.50 


Decision Point 2 Time (seconds): 


Mean= 2.87 Standard Deviation= . 08 
Minimum= 2.74 Maximum= 3.00 
Median= 2.86 
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Table IV-3 
The statistics of the following parameter are the results of twenty 
production runs from the same initial conditions. 
Time to go from decision point 1(seconds): 
Mean= 1. 56 Standard Deviation= .18 
Minimum= 1. 30 Maximum= 1. 82 


Median= 1.49 


Time to go from decision point 2(seconds): 


Mean= « 25 Standard Deviation= .08 
Minimum= .12 Maximum= a8 
Median= . 26 
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Table IV-4 
Termination Time Time To Go Time To Go 
From Point 1 From Point 2 
(seconds) (seconds) (seconds) 
3,12 1.76 38 
3.12 1.57 . 295 
3.12 3 ra 29 
3.12 1.68 14 
3.12 1.79 .13 
3.12 1.81 . 36 
3.11 1.40 13 
3.11 1.32 «29 
3.12 1.82 . 24 
3.11 1.35 .16 
3.11 1.47 14 
3.12 1.75 20 
3.11 1.30 .23 
3.12 1.74 12 
3.11 1.37 Pe 
3.11 1.47 . 29 
3.11 1.40 . 33 
3.11 1.48 . 26 
3.12 1.50 27 
3.11 1.41 35 
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VI. Conclusions and Recommendations 


Conclusions 

The objectives of the simulation were listed earlier. All 
objectives were accomplished to varying degrees. Subroutine PRINT 
produced the results of interest. 

The first objective was good aircraft performance. The simu- 
lated performance of the aircraft in the selected maneuvers compares 
well with the actual aircraft performance as determined from the 
author's flying experience and the F-4C Technical Order. This was a 


major objective of the simulation. 


A second objective was to incorporate pilot judgement into the 
simulation. As the printed results indicate, the selected maneuvers 
were initiated at various ranges producing the simulated judgement 
error of the human pilot. The results indicate that the range at which 
the pilot executes a maneuver has a strong effect on the success of the 


missile attack. 


The program has the ability to select the desired maneuver at 


predetermined decision points in the air-to-air combat simulation. 


¥ The new maneuvers are selected at the ranges determined by the pilot 
judgement subroutine. The inputs for these maneuvers were updated 
during each integration interval. The inputs are equivalent to the pilot 
control] inputs to the ailerons, elevators, and rudders. In addition, 
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the aircraft has the ability to maneuver or fly straight and level 


Therefore, the third objective was met. 


The aircraft model has provisions to incorporate most maneu- 
vers applicable to air-to-air combat. The maneuvers selected for the 
program are taken from the latest tactical flying manuals and informa- 


tion from the staff of the TFWC. 


Recommendations 

Further development of complete aerodynamic coefficient 
equations will more completely simulate the performance of the air- 
craft. An aeronautical engineer with flying experience in fighter type 
aircraft should develop a more complete model to produce an authentic 
maneuvering target. 

The maneuvers must be continually changed and updated so that 
a true test can be made of missile guidance. The latest tactical 
maz::uals and the staff of TFWC are the best sources to update the pro- 
gram maneuvers. 

All the missile routines must be developed for complete combat 
simulation. The routines are available in the present program but in- 
clude only a very basic missile model used just to test the target simu- 
lation. 

The development of the missile model and missile noise is 


beyond the scope of this thesis. The missile seeker noise programs 
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are elementary attempts at modeling noise as a function of relative 


position. Further study in this area may provide a more realistic 


simulation. 
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Appendix C 


/ Aircraft Aerodynamic Coefficients and Equations 
Equations 
Lift: (37) 
C,; = CEF 1+ CEF 2(4a)+ CEF 18 (A MN) + CEF 13 ( 6 s) 
Drag: (38) 
Ch = CEF 3+ CEF 4 (Aq) + CEF 19 (4 MN) 
Sideforce: (39) 
Cy = CEF 23 (48) + CEF 17 (48) 
Pitching Moment: (40) 
ee = CEF 5 + CEF 6 (4a) + CEF 7 (Aq) (3) (Cr + CEF 8 (Aq@) 
(2) (=m) + CEF 20 (4 MN) + CEF 14 (5s) 
Roliing Moment: (41) 
Cy, = CEF 21 (48) + CEF 9 (A p) (2) a + CEF 10 (A r) (3) 
(a) + CEF 15 (48 ) 
Yawing Moment: (42) 
Cy = CEF 22 (Af) + CEF 11 (dr) (3) Gr) + CEF 12 (4 p) (4) 
Ge) + CEF 16 (48) 
Coefficients are listed in Table III. 
Symbols: 


Sa - delta alpha (angle of attack) - degrees 
Op - delta beta (sideslip angle) - degrees 


b - wing span - feet 


: us " Iie sate cisammacememensimmnaiasiininNass 


Ar - 


ie ROL ay te ite a AO 


mean aerodynamic cord - feet 

delta alpha dot (time rate of change alpha) 
delta p (rate of roll) 

delta r (rate of yaw) 


delta q (rate of pitch) 


—_——— eae 


Table V 


Stability Derivatives and Coefficients 


Program Standard 
Symbol Notation Definition 
CEF 1 Cro Basic Lift Coefficient 
CEF 2 Lee Derivative of L respect to a 
CEF 3 Cho Basic Drag Coefficient 
CEF 4 Crg Derivative of D respect to a 
CEF 5 Co Basic Pitching Moment Coefficient 
CEF 6 Cc ma Derivative of m respect to o 
CEF 7 Cu Derivative of m respect to q 
CEF 8 Cis Derivative of m respect to @ 
CEF 9 Cip Derivative of 1 respect to p 
CEF 10 Cie Derivative of 1 respect to r 
CEF 11 Cnr Derivative of N respect to r 
CEF 12 Cnp Derivative of N respect to p 
CEF 13 Cribs Derivative of L respect to 6s 
CEF 14 Cmés Derivative of m respect to 6s 
CEF 15 Cir Derivative of 1 respect to 6r 
: CEF 16 Cnor Derivative of N respect to 6r 
; } CEF 17 Cyo6r Derivative of Y respect to 6r 
CEF 18 Cim Derivative of L respect to M 
¢ CEF 19 Com Derivative of D respect tom 
| 164 
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: Program 
Symbol 


CEF 20 
CEF 21 
CEF 22 


CEF 23 


Key Symbols 
L - Lift D - 


a - angle of attack 


@ - rate of change of a 


p - rate of roll 
r - rate of yaw 


1 - rolling moment 


6s - elevator deflection 


q - rate of pitch 
m - pitch moment 


N - yawing moment 


§r - rudder deflection 


Standard 
Notation Definition 
Cum Derivative of M respect to m 
Cis Derivative of 1 respect to B 
Cnp Derivative of N respect to A 
Cc ye Derivative of Y respect to 8 
Drag M - Mach Number 

6 - sideslip 


Appendix D 


Missile Model 


To test the performance of the aircraft simulation, a missile 
model was required. The model included in the program was a very 
simplified and generalized missile. The missile was a proportional- 
navigated-guided model with constant velocity, constant mass, and it 
was unaffected by gravity. Missile specifications were generalized 
from state-of-the-art air-to-air missiles. The specifications of the 
missile were: 

1. Proportional navigation constant - 3 

2. Time constant - .5 seconds 

3. Gimbals limits - 30 degrees 

4. "G' force - 20 g's. 

5. Line-of-sight rate - 30 degrees/second. 

During actual missile simulations using the aircraft model, 
the missile model should be developed along the same modular con- 
struction as the aircraft model. All the required subroutines are 
present in the Program EVASION but are very basic due to the simpli- 
fied missile. 

The missile model was constructed so that any type guidance 
could be simulated. If other than infra-red or radar guided are 


desired, noise source subroutines should be developed. 
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ee eee — 


( Additional complete subroutines must be developed for an 

| accurate missile model. Differential equations are required for the 
equations of motion routine. The missile model can be made as com- 
plete as desired by the complexity of the total coefficient and force 


equations in FORCESM. Additional parameters can be put into the 


model in DATAMIS. 


FORCESM 
In the development of the missile model, the total velocity of 


the missile is assumed to be along the longitudinal axis, X-axis. 


<| 
~* 


(43) 


This assumption in turn leads to the angular velocity, w, being 


equal to vector sum of the rate of change in the azimuth and elevation. 


O O 
® # Ww =| ¥ (44) 
y . 
j o,. ¥ , 
& Lens 


The missile is assumed to be of constant velocity. The accelera- 


tion in the aircraft navigation is therefore equal to the cross product of 


t ‘ 
1 ; the angular velocity and velocity in the body fixed frame. 
‘ 
Se: - Go} Vy 
Vn Oe « v,.* N x |O (45) 
| bi O 
fig 167 


a bjs o 4 
Ye = |o xe = ¥ Vj (46) 
Ry OO =A Vis 


(i, j, k) - unit vectors 


To determine the acceleration forces in terms of g's, the lateral 


and vertical acceleration are divided by the gravity. 
os Mo 2 ee 
oe (a7 
g 


A test is needed in FORCESM so that the missile will not exceed 
its maximum g forces. If the computed g forces are greater than 
specif ications, the rate of change of azimuth and/or elevation will be 
recomputed to bring the g forces within limits. INPUTSM is called 


in FORCESM to recompute the azimuth and elevation angles. 


168 


Appendix E 


Entering Data 


Aircraft and missile data are stored in arrays for ease of pro- 
gramming. The data must be read into the program ina specific 
order for proper operation. The data for the aircraft is put in a file 
called TAPEI and the missile's data is on TAPE 2. In this manner, 
selection of one vehicle is not dependent on the other. 

After the initial integer which indicates the type missile and 
target, the vehicle data is stored in arrays AC (aircraft) and DM 
(missile). 

The data is entered one value per card in free format in the 
following order for the aircraft: 

LP Weight 

2 Wing Area 

3. Wing Span 


4. Mean Chord 


5. Maximum Rate of Change of Thrust 

6. Maximum Rate of Change of Bank 

Na Maximum Angle of Attack 

8. Maximum Sideslip Angle 

9. Maximum Rate of Change of Angle of Attack 


10. Maximum Rate of Change of Sideslip 


ee ——— = 


= 


ET. Product of Inertia XX 
E22. Product of Inertia YY 
ine Product of Inertia ZZ 


14. Product of Inertia XZ 


The next information to be entered is the initial conditions 
which are stored in arrays DATAFT (aircraft) and DATAMIS (missile). 
The initial conditions are entered one value to a card in free format 


in the following sequence: 


Lie x - position ie bank 

Ze y - position 8. pitch 

3. z - position os heading 
4. u - velocity 10. yaw rate 
5. v - velocity it. pitch rate 
6. w - velocity £2. roll rate 


The thrust coefficients are entered in free format, three 
values per card, as a function of Mach number while the aerodynamic 
coefficients are entered, six values to a card in free format as a 
function of angle of attack as well as Mach number. Interpolation 
is used by the table look-up routine to select the proper coefficients 
for the given parameters. 

The missile model and its parameters will determine the 


order of entry of its data. The arrays are specified in the program. 
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